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APERTURE EXCITED DIELECTRIC ANTENNAS 

By W. F. Croswell, J. S. Chatterjee,* V. B. Mason,** and C-T. Tai** 

Langley Research Center 

SUMMARY 

The results of a comprehensive experimental and theoretical study of the effect of 
placing dielectric objects over the aperture of waveguide antennas are presented. Experi- 
mental measurements of the radiation patterns, gain, impedance, near-field amplitude, and 
pattern and impedance coupling between pairs of antennas are given for various Plexiglas 
shapes, including the sphere and the cube, excited by rectangular, circular, and square 
waveguide feed apertures. The waveguide excitation of a dielectric sphere is modeled 
using the Huygens’ source, and expressions for the resulting electric fields, directivity, 
and efficiency are derived. Calculations using this model show good overall agreement 
with experimental patterns and directivity measurements. The waveguide under an infinite 
dielectric slab is used as an impedance model. Calculations using this model agree quali- 
tatively with the measured impedance data. 

It is concluded that dielectric loaded antennas such as the waveguide excited sphere, 
cube, or sphere- cylinder can produce directivities in excess of that obtained by a uniformly 
illuminated aperture of the same cross section, particularly for dielectric objects with 
dimensions of 2 wavelengths or less. It is also shown that for certain configurations 
coupling between two antennas of this type is less than that for the same antennas without 
dielectric loading. 


INTRODUCTION 

Many of the future communications or navigation satellite systems will require 
multiple or adaptive beam antennas with shaped or controlled beam coverages, depending 
upon the application. Common approaches to solving this problem include the placement 
of multiple feeds in the focal plane of a direct-fed parabolic reflector, or to assemble 
the same feeds to illiuninate an offset parabolic reflector. In the case of the direct-fed 
parabolic reflector the size of the multiple beam cluster will cause blockage, and for both 
the direct-fed and offset-fed parabolic reflectors the spillover and side lobe level of each 
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beam is dependent upon the amount of interelement coupling and how closely packed 
the feed elements can be placed. Also important for this type of design are minimum 
spillover consistent with achieving high aperture efficiency, which is again related to 
the size and packing capability of the multiple beam feed elements. As a result, there is 
a requirement for a feed element that exhibits low blockage, low mutual coupling between 
adjacent elements, and that also produces the maximum directivity for a given ^erture 
size. The work described in this report was motivated by this requirement. 

It has been found experimentally that an antenna constructed by placing small dielec- 
tric objects such as the sphere or cube over a waveguide aperture exhibit radiation and 
coupling properties desirable for feed antennas. In this report many of these measure- 
ments are presented. Subsequently the waveguide aperture excitation is modeled using a 
single Huygens' source, and a rigorous electromagnetic analysis is presented for the 
dielectric sphere. These theoretical results are found to agree closely with the measure- 
ments. Finally , the results of the complete study are discussed and concluding remarks 
about the aperture excited dielectric (AED) antenna and its properties are presented. 

SYMBOLS 


a 


b 


b 


1 


C 


e 


C 


m 


D 


d 


D 


0 


ii -plane 


E(R) 


sphere radius 

radius of circular waveguide 
normalized susceptance 

distance of Huygens' source to center of coordinates 

electric current moment 

magnetic current moment 

sphere diameter, 2a 

thickness of dielectric slab 

directivity 

pattern measurement plane which contains electric-field vector 
electric field at position vector R 
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E 


electric field 


f 

g 

Gq(R|8') 
g{nm) (r| r*) 

^^nm) (p| r.) 
H -plane 
H(R) 

H 

hj^OcR) 

& 

I 

Iq 

j(R’) 

Jl(x) 

i„(kR) 

-i 

M 


frequency 

normalized conductance 
free-space dyadic Green's function 
dyadic Green's function of third kind 
dyadic Green's function of fourth kind 

pattern measurement plane which contains magnetic-field vector 
magnetic field at position vector R 
magnetic field 

spherical Hankel function of first kind 
unit dyad 

constant current magnitude 

electric current density 

cylindrical Bessel function of first kind 

magnetic current density 

spherical Bessel function of order n 

free-space wave number, 277/ \ 

wave nvunber of medium n, 

length of electric current s^ment 

spherical vector wave functions 



N 

N 


n 


(cos d) 
R 


R' 


Rr 

tan S 


TE 


01 


TE 


11 


W 


xl 


11 


X, y, z 


X,Y,Z 


Y. = Y + V 
^in ^TE ^ ^TM 


'TE 


TM 


spherical vector wave function 
index of refraction, 

unit normal directed away from body surface 
dissipated power 

associated Legendre function of degree m and order n 
vector from coordinate origin to observation point 
vector from coordinate origin to source point 
radiation resistance 
loss tangent 

E mode in a rectangular or square waveguide 
E mode in a circular waveguide 
total radiated power 
first root of — [ji(xji)J 

coordinate directions 

near-field-measurement coordinate directions 

normalized aperture admittance, TE vector potential 
normalized aperture admittance, TM vector potential 
normalized radial wave number 
reflection coefficient 
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1 (m = 0) 



(m / 0) 


S{R-R*) three-dimensional delta function 

fl {n=-t) 

to 


f 

€ 

real part of dielectric constant 

TT 

e 

imaginary part of dielectric constant 

e 

dielectric constant of medium n 

n 



relative dielectric constant 

'n 

n 

wave impedance of medium n 

^0 

free-space wave impedance, ^/mqAo 

6 

polar angle 

\ 

free-space wavelength 


relative permeability constant 

(j 

conductivity 

<P 

azimuthal angle 

CO = 2Trf 



ABBREVIATIONS 

AED 

aperture excited dielectric 

VSWR 

voltage standing wave ratio 

VTVM 

vacuum tube volt meter 
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EXPERIMENTAL MEASUREMENTS 


Radiation Patterns for Plexiglas Spheres and Cubes 


The aperture excited dielectric (AED) antennas are essentially three-dimensional 
sources formed by a complex field distribution over the surface of the feed waveguide and 
the dielectric surface. The radiation patterns therefore are sensitive to the type of feed 
waveguide; size, shape, and dielectric constant of the dielectric object loading the aperture; 
operating frequency; and the relative spacing between the feed waveguide and the object. 

Due to its availability and low cost, all dielectric antennas used in the experimental meas- 
urement reported in this paper were constructed from Plexiglas (ej. = 2.57; tan 8 = 0.0065). 
As shown in later sections using the Huygens’ source near a sphere as a model, other 
dielectrics should also produce directive antenna patterns. 

The results of pattern measurements of a variety of waveguide types exciting different 
dielectric objects are described in this section. The following feeds, detailed dimensions 
of which are given in figure 1, were used: 

(a) Circular waveguide, C-band, excited in the TEj^j mode 

(b) Slot at the end of a rectangular waveguide, X-band, excited in the TEqj mode 

(c) Rectangular waveguide, X-band, excited in the TEqj mode 

(d) Square waveguide, X-band, excited from a TEqj mode rectangular waveguide 
with an E- plane taper 

A photograph of the circular waveguide feed with Plexiglas objects is shown in figure 2. 

The radiation patterns were measured using these feeds with Plexiglas cubes and spheres 
of various dimensions over wide frequency ranges. Only typical patterns illustrating the 
important features of AED antennas are given in succeeding paragraphs. 

Circular waveguide feed .- The open-end circular waveguide shown in figure 1(a) was 
used to excite a 5.08-cm Plexiglas cube and 5.08- , 7.62-, 10.16- , 12.70-, and 15.24-cm 
Plexiglas spheres over a frequency range from 5.0 to 6.0 GHz. In each instance the dielec- 
tric obstacle was placed tightly against the waveguide aperture, which in the case of dielec- 
tric spheres resulted in part of the sphere protruding inside the feed waveguide. Patterns 
measured at 5.0 GHz for selected dielectric objects are given in figures 3 and 4 for E- and 
H-plane cuts, respectively; open-end waveguide patterns are given for comparison. It was 
noted that the cross-polarized pattern was about the same level as that for the open-end 
waveguide and was therefore not included in the AED radiation pattern plots. 


Slot feed .- A thin slot excited by a rectangular waveguide (fig. 1(b)) was used to 
excite a 5.08-cm Plexiglas cube and Plexiglas spheres 5.08, 7.62, 10.16, 12,70, and 15,24 cxr 
in diameter. A special- shaped Plexiglas object designed to produce a uniform phase front 


as determined from ray optics was also used. Resulting patterns measured at 10.0 GHz 
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are given for selected cases in figures 5 and 6. Also included in figure 5 are the patterns 
of the open slot without a dielectric obstacle. It is noted that the E-plane patterns for the 
slot AED antenna are poor for the smaller dielectric objects, while the H-plane patterns 
for dielectric objects of identical size display improved directivity. 

Rectangular and square waveguide.feeds .- Radiation patterns were measured at 
,10.0 GHz for the 5.08-cm Plexiglas cube and the 5.08-, 7.62-, 10.16-, and 12.7-cm spheres 
using the rectangular waveguide and square-horn antenna shown in figures 1(c) and 1(d). 
Patterns for selected cases are given in figures 7 to 10. As expected, the H-plane patterns 
for the two feeds were very similar. The E-plane patterns were also similar over the main 
lobe down to about 10 dB, with the horn patterns exhibiting generally lower side or back 
lobe levels. Overall, the horn patterns are cleaner; for example, they show little evidence 
of phase error. 


Influence of Feed Displacement 

The effect of axial feed displacement on the radiation pattern was studied using the 
circular waveguide antenna operating at 5.0 GHz with a 12.70-cm Plexiglas sphere. The 
sphere was displaced a distance 0.635 cm (0.106\) and 1.27 cm (0.212\) from the plane of 
the circular waveguide feed. The resulting radiation patterns are given in figure 11. The 
primary effect of the displacement was to lower the directivity by broadening the E-plane 
radiation patterns and by increasing the side lobe energy along the aperture plane. Meas- 
urements of the effects of displacement were made using other Plexiglas obstacles, with 
similar results being obtained. 

Effects of Object Dimensions on Radiation Patterns 

In general, it was noted that the radiation patterns for spheres and cubes became 
more directive with increased object size. The effects of the dimensions of the Plexiglas 
AED antennas on the radiation patterns for other shapes were investigated usir^ Plexiglas 
rectangular blocks and cylinders with spherical tips. This study was made using the cir- 
cular waveguide antenna operating at 5.0 GHz. Photographs of the cylindrical and rectan- 
gular blocks are given in figure 12. The radiation patterns made using rectangular blocks 
with a base dimension of 5.08 cm square and with heights varying from 2.54 to 20.32 cm 
are given in figure 13. The cylindrical AED antennas were constructed by bonding together 
5.08-cm cylinders of various lengths with 5. 08 -cm -diameter hemispheres. Patterns of this 
AED antenna with cylinder lengths of 0.0, 1.27, and 2.54 cm are given in figure 14. It is 
observed from figure 13 that there appears to be a single rectangular block configuration 
(i.e., the cube) that gives good directive radiation patterns with low side lobes. From fig- 
ure 14 It can be observed that the hemispherically tipped cylinder produces symmetrical 
radiation patterns with the directivity increasing as a function of length; also, the E-plane 
and H-plane patterns are nearly identical. 
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The top surface of a 12.7-cm-diameter sphere was shaped in a manner to produce 
a uniform path length as explained in figure 15(a); a photograph of the spherical and shaped 
object is given in figure 15(b). The circle APMONQA represents the surface of the sphere. 
The dotted line PMONQ represents the portion of the sphere which has been removed 
to form the surface represented by line PLORQ. The shaping is such that any typical ray 
path ABC from the source A reaches the front surface SOS’ in the same interval of 
time, that is, in phase. The patterns of the circular waveguide antenna exciting this object 
in the direct and inverted positions in the waveguide aperture are given in figure 16. Upon 
comparing the patterns given in this figure with earlier ones using spheres, it is concluded 
that for these electrically small dielectric obstacles the shape of the outer surface away 
from the aperture has little effect upon the main lobe structure of the radiation pattern. 

For the AED antennas with the dielectric object placed directly upon the aperture of 
the feed waveguide, all of the shapes produced directive patterns which in most cases had 
low side lobes. Particularly for the objects where the dielectric dimensions were afree- 
space wavelength or less (such as the cube, hemispherically tipped cylinder, or small 
spheres) the E- and H- plane patterns were nearly identical over the main lobe. 

Beamwidth and Gain 

In order to indicate how the patterns of the AED antennas vary as a function of dielec- 
tric dimension, the 3-dB beamwidth is plotted against normalized sphere diameter for the 
spheres excited by various waveguide feeds, as shown in figure 17. Also shown in this 
figure is a straight line corresponding to the beamwidth of a uniform circular aperture. 

It is noted that all spheres exhibit beamwidths narrower than the uniformly excited aperture 
of the same diameter. 

The gain of the circular waveguide antenna of figure 1(a), exciting various Plexiglas 
objects, was measured at 5.4, 5.5, and 5.6 GHz; the measured values are given in table 1. 
Any input VSWR for each dielectric object at each frequency was accounted for in the 
measurements. In figure 18, the gain of several spherical AED antennas is compared 
with that of an optimum horn havii^ the same cross section. It should be noted that the 
gain of the AED antenna is at least 3 to 4 dB higher than the corresponding horn. 

Impedance 

The impedance of the AED antennas in general depends upon the same parameters as 
the patterns, for example, the type of waveguide feed, size and shape of dielectric object 
loading the aperture, frequency, etc. However, in previous work (ref. 1), it was determined 
experimentally that the impedance of a 1 to 2 wavelength slab covering waveguide fed aper- 
ture antennas with flanges about the same size as those used in this study was about the 
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same as that calculated for the waveguide opening onto an infinite slab -covered ground 
plane. Since there is no available theory to treat the geometry of the AED antenna, and 
at the present time little hope of attaining one, this model will be used as a first 
approximation. 

Comparison with available theory .- The admittance of a circular waveguide opening 
onto an infinite ground-plane covered by a semi-infinite dielectric slab has been derived 
by Bailey and Swift (ref. 1). Consider the geometry given in figure 19. The expressions 
for the input admittance of the circular waveguide excited by the TEj^j mode, in this 
case, are 


Yin = 


TE 


+ Y, 


TM 


( 1 ) 
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where is the radius of the circular waveguide, /S is the normalized radial wave 
number, N is the square root of the complex dielectric constant, is the first root 

AfjJxlJ, d is the thickness of the dielectric slab, and is the wave number 
“ dx ^ ^ 

in free space, 27r/A.. These impedance expressions were programed, and calculations 
were performed for the waveguide dimensions given in figure 1(a) with dielectric slabs 
having the same thickness as the diameter of the sphere or side dimensions of the cube. 
These calculations were made over a frequency range of 5.0 to 6.0 GHz and compared 
with admittance measurements for the cube (see fig. 20). It should be noted that the aper- 
ture admittance was measured at spot frequencies for the spheres and the data showed a 
general admittance which was approximately real and lower than the admittance of the 
cube. Measurements were also made of the input VSWR of the same aperture loaded with 
spheres; the measurements are compared with the calculations in figure 21. The agree- 
ment for all measurements is quite good considering the crudeness of the model. There- 
fore, the slab model should, for practical purposes, be quite useful for predicting the input 
admittance of AED antennas, regardless of the shape. 

Matching devices .- Because of the smooth behavior of the impedance and VSWR curve 
measured for AED antennas as a function of frequency, it was believed that plugs made of 
similar material about a wavelength or so lor^ would be good impedance matching devices. 
Three plugs with various straight and tapered sections were designed with the dimensions 
and shapes described in figure 22. The VSWR measured using the plugs inserted in the 
waveguide are summarized in table 2. Similar results were obtained with larger spheres. 

It should be noted that these measurements were made with a coaxial slotted trans- 
mission line and hence required a coax to waveguide junction. This junction was optimized 
for 5,4 to 5.6 GHz, As a result of this matching study, it is believed that AED antennas of 
other sizes and shapes may also be matched if tapered plugs are used. 


Near- Field Measurements 


In order to gain some insight into why the AED antennas gave the improved directivit 
indicated by the pattern and gain measurements, it was decided to measure the near-field 
properties of several antennas. The simplest method, which measures amplitude only, 
was devised by Rudge (ref. 2) to measure leakage through small cracks in microwave oven 
doors. The measurement probe consists of a microwave diode whose pigtail leads are 
used as a probe dipole. Attached to the dipole is a detection circuit as shown in figure 23, 
The wings of the dipole are trimmed until, as Rudge demonstrates, the coupling between 


the aiiLenna to be measured and the dipole is small. lu tile case oi ArjiJ aillennas 


. ^ ^^2. xVii-k 1 <-\vi riT 

V r. J J itllLKllIliia LiiC 


of the dipole wings was trimmed until the dipole could be inserted into the aperture (and 
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partially into the waveguide itself) without any change in the VSWR observed at the input 
terminals of the waveguide. The dipole probe was mounted in a thin glass tube attached 
to a three- axis precision lathe head, placed some distance below the antenna probe plane. 
The antenna and probe were then placed inside a box lined with microwave absorber to 
minimize reflections. 

The probe measurements were made in a plane 0.508 cm in front of each antenna 
being measimed. Although complete probe measurements were made in fine increments 
over the entire aperture of all antennas measured, only principal plane cuts and one 
measurement down the side of each antenna are presented. The geometry of the measure- 
ments is given in figure 24. The near-field amplitude measurements for the different 
Plexiglas dielectric objects fed by a circular waveguide antenna operating at 6.0 GHz are 
given in figures 25 to 31. It is readily apparent that near field of significant amplitude 
extends considerably beyond the physical limits of the dielectric object. The fields in 
general are tapered and therefore may, in part, account for the low side lobes in the far 
radiation patterns of the same AED antennas. In the case of the spheres, the measured 
near fields along the sides of the dielectric objects (denoted as the Y-direction in fig. 24) 
diminish rapidly as the distance equal to the lower hemisphere is reached. This indicates 
that for many shapes the field is concentrated primarily in the forward part of the dielectric 
object, away from the feed aperture. This field concentration indicates there should be 
reduced coupling between adjacent antennas of this type, as compared with open-end wave- 
guides. This decoupling is confirmed by the results presented in the next section. 

Coupling 

Pattern coupling .- In order to determine how the patterns of an AED antenna are 
modified by the presence of other elements, radiation patterns were made using a pair of 
circular waveguide antennas of the type shown in figure 1(a) where one antenna was excited 
and the other antenna was terminated in a matched load. The antennas were located side 
by side with the excitation probes parallel so that either E-plane or H-plane excitation 
could be obtained by rotating the feed waveguide by 90*^. 

Pattern measurements were made for pairs of AED antennas, including the 5.08-, 
7.62-, 10.16-, and 15.24-cm spheres, as a function of the spacing between elements. Spac- 
ings were chosen so that at the minimum distance the waveguide aperture or Plexiglas 
sphere touched. Measured results are given in figures 32 and 33 for the open-end wave- 
guide and 5.08-cm sphere, respectively. For the open-end waveguide, the E-plane patterns 
are modified more than the H-plane patterns by the presence of the parasitic element. For 
the sphere-loaded apertures, however, the H-plane appears to have more pattern coupling 
than the E -plane. Similar results were obtained for the other dielectric objects. 
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Impedance coupling .- Impedance coupling measurements were made with two 
waveguide configurations. First, coupling was measured as a function of spacing with the 
circular waveguide configuration given in figure 1(a). Similar measurements were made 
using the same size waveguides mounted in a continuous ground plane 30.48 by 60.96 cm. 

A photograph of the ground plane and antennas is given in figure 34. Measurements taken 
of the circular waveguides with flanges are given in figure 35 for the open-end waveguide, 
the 5.08-cm cube, and the 5.08-cm and 7.62-cm spheres at an operating frequency of 
5.8 GHz. It is noted that the E-plane coupling is much lower than the H-plane coupling, 
which is not the usual case. The coupling measurements of the circular waveguides opening 
onto a continuous ground plane are given in table 3. These measvurements include both 
E- and H-plane coupling of the two waveguides arranged as shown in figure 35. 

From the data in table 3, it is noted that for the continuous ground plane measure- 
ments the E-plane coupling is higher than the H-plane coupling for the unloaded case. 

This is consistent with other ground plane mounted waveguide calculations and measure- 
ments. However, with the addition of the dielectric objects, the coupling in the E-plane 
drops markedly to about the same as the coupling in the H-plane. 

ELECTROMAGNETIC THEORY OF A DIELECTRIC SPHERE 
EXCITED BY A HUYGENS’ SOURCE 

Huygens' Source Model 

In order to understand and optimize the AED antenna, it is desirable to obtain a good 
theoretical model for the impedance and radiation pattern of the waveguide radiating in the 
presence of a dielectric obstacle. The solution of this problem for the detailed feed and 
dielectric geometry is exceedingly difficult. It is well known that the Huygens' source, 
the superposition of an electric and magnetic current source, can produce a cardioid patterr 
when radiating into free space. This cardioid pattern closely resembles the radiation 
pattern of waveguide antennas having small ground planes. Also, the simplest three- 
dimensional dielectric obstacle for which Green's functions are available is the sphere. 
Therefore, the Huygens' source model composed of the superposition of electric and 
magnetic dipole cvirrents near or in a dielectric sphere was chosen as a good theoretical 
model. A more detailed discussion is contained in reference 3. 

Figure 6 shows the geometry under consideration. The source consists of an 
x-directed electric dipole, together with a -y-directed magnetic dipole. The electric 
and magnetic current moments are denoted by Cg a-nd Cm, respectively. The regions 
exterior and interior to the sphere are characterized by different constitutive param- 
eters and, therefore, different propagation constants kj and k 2 , where 
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^1 ■ 


(Region I; R > a) 


(4) 




" \/^2 ^2 


(Region II; R < a) (5) 


Here, ^2 ^2 are the permeability and permittivity of the sphere, which may be 

complex. For this application, region I will be free space, and hence = Mq and 

" ^ 0 * 

With C = VnC (where Vn is the free- space wave impedance) the currents 
m u e ' u 

oriented as shown in figure 36, in the absence of the sphere, form a Huygens’ source. It 
radiates a linearly polarized, cardioid pattern with a field maximum in the -z -direction 
and a null in the +z-direction. It is because of this directional characteristic, which is not 
unlike that of an open-end waveguide in free space, that the Huygens' source was chosen 
to model the waveguide excitation of the dielectric sphere. The source can only be approxi- 
mated in the presence of the dielectric sphere because the complex wave impedance -q is 
not known. 

Although, in general, v would be a function of sphere diameter and constitution as 
well as source location,- the following current relations are used to construct the source 

I 

' = I (b,sa) (7) 

where Vi = and -r? g = Equation (6) was chosen because it would be valid 

if bj >> a, since for this case the sphere would have little effect upon the source. Alter- 
nately, equation (7) would be appropriate for the source located within or on the surface of 
a large lossy dielectric sphere. Also, it has been found that equation (7) yields a good 
approximation to the Huygens’ source on the surface of small spheres. 

Dyadic Green's Functions Pertaining to the Dielectric Sphere 

Under the assumption that the current distribution is a known function of position, the 
radiation from a dielectric sphere may be investigated as a boundary- value problem and 
the solution may be obtained by the technique of dyadic Green's functions. 

Figure 37 shows the geometry of a dielectric sphere with an electric current source 
located at R' • With the source located outside the sphere in region I, the electric fields 
in regions I and II may be expressed by 
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E{R) = Wj j jj (RiR’) • J(R’) dV’ 


E{R) = jjj (R| R’) . J(R') dV' 


Similarly, with the source located inside the sphere 


E{R) = io>^ JJJ 

E(R) = IJJ 


(R[ R') ■ J(R') dV’ 


(R| R') ♦ J(R') dV’ 


(Region I) (8) 


(Region II) (9) 


(Region I) (10) 


(Region n) (11) 


is the dyadic Green’s function of the third kind and the harmonic time 

been suppressed. The superscripts m and n of the Green’s functions 
in equations (8) to (11) denote the regions of the observation and source points, respectively 
Those functions are solutions to the vector wave equations 


where Gg(R|R’) 
factor has 


V XV X G3-k2G3=lS(R-R’) 

vx VXG3-k^G3=0 

Besides satisfying the radiation condition at infinity, 

R “ ^ ^3 I - ikR X G3(R I R’)J = 0 

G3 also satisfies the boundary conditions at the surface of the sphere; that is, 
nxG3(R!R/)L^^^ . 

1 ‘ = -I- 1 1 

_nXV X G 3 (r|r’) =_ nX vxG«(RjR')| 

1 >R=a_ M2 3 I ^1 

+ j 


(Ri a) 
(R s a) 


(12) 


(13) 
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The simplest way to find G 3 (R| R’) is to use the method of scattering superposition 
in which the fields are assumed to consist of a sum of incident and scattered waves. With 
this method the form of the scattered waves are first constructed using eigenfunctions 
which will satisfy the boundary conditions at the sphere surface, as well as having the 
proper form at R = 0 and R = ® . The coefficients of the scattered waves are then 
determined by matching the boundary conditions on the surface of the sphere. 

Thus § 3 (r| R’) is treated as consisting of two parts: 

(r|r’) = Gq(r1r*) + GggCRiR') 

G^21) (R| Rt) = (R| R’) (R 5 a) (15) 


(R i a) (14) 


Here, Gq(R| R') denotes the free-space dyadic Green’s function p^rtainii^ to an infinite 
region with the same constitutive parameters as region I, Gjg (Rl R') represents the 

part of the wave scattered from the sphere, and G^g^^ (R] R') is used to denote the Green's 
function internal to the sphere since only a scattered field will be present in this region. 
The free-space Green's function for this problem is given by Tai (ref. 4) as follows. 




n 


n=l m=0 


(ki)N< (k,) 


,mn 


mn 


2n + 1 (n - m)! 


n(n + 1) (n + m)! 


(k^ m; (kj) 


^mn 


mn 


(R > R') 


> (16) 


^ fcl ifeo n(n + l)(n + m)! 


M. 


:mn 


(kj) (k^) 


mn 




;mn 


;mn 


(R < R') 
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where 




(m = 0) 

0 ) 


and 


M_ (k) = j (kR) 

Qmn 


m 


^sin 


sin^ 


^)cos - 


BP“ (cos e) 
n * cos 


^0 


sin 




(17) 


N. (k) = ililLLl) j^(kR) p“ (cos 9) 


;mn 


kR 


n 


sm 


mc?5>R 




^pIT(cos e) COS ^ 

xxicpO "I 


do 


sin 


sin (9 


(18) 


The two sets of spherical vector wave functions represented by M and N are 
solutions to the homogeneous vector wave equation 


V X V X F - k F = 0 


as shown by Stratton (ref. 5), and also satisfy the symmetrical equations 


N. (k) =-L V X (k) 


;mn 


0 


mn 


(19) 


M (k) V X N (k) 
ginn k gmn 


(20) 


In equations (17) and (18), jjj(kR) denotes the spherical Bessel function of order n 
Pjj (cos 6) represents the associated Legendre function of degree n and order m. 
The expressions in equation (16) with superscript (1) are obtained from equations (17) 
and (18) by replacing^the spherical Bessel functions with spherical Hankel functions of the 
first kind, that is, h^-' (kR), A prime used on a function denotes that it is defined with 
respect to the primed coordinate variables {R',9\<p') pertaining to the vector R'. 
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In constructing G 3 (R| R’) from the wave functions, the following observations are 

made: 

(1) For the source located exterior to the sphere, the posterior parts of G 3 y"^(R| R') 
and G^g^^ (r|r*) must be the same as Gq(R| R') evaluated at R < R'. This 
is necessary to match the boundary conditions at the sphere surface. 

(2) The propagation constants kj and ^ must be used in the anterior parts of 

R') and (R| R')> respectively, because of the regions in which 

the observation points are located. 

(3) The spherical Hankel functions must be used in the anterior parts of G 3 g^^(R|R') 
because these functions represent outward traveling waves. 

(4) Finally, since the fields are finite at the origin, the anterior parts of Ggg^^(RlR’) 
must be constructed using the spherical Bessel fxmctions. 

In keeping with these requirements, let 



(RiR .) .hf r (2 - V 

^ ^ n(n + l) (n + m)! 


aj.^) (ki)M;(^^ (kj) 

1^ Qmn Qmn Qmn 


b(l) (kj) Os.J 


Qmn Qmn 


mn 


(21) 


and 


E-v... 


00 n 


2n 


+ 1 (n - m) ! r (1) 
+ 1) (n + m) ! ®m: 


“e„ “el’ 


mn Qmn 


mn 


Qmn Qmn 


(ko) 


Qmn J 


( 22 ) 


The scattering coefficients a, b, c, and d in equations (21) and (22) are complex 
quantities and represent the magnitude and phase of the sphere’s contribution to the total 
fields. The superscripts (1) on these coefficients denote that they are defined for the 
source located in region I (fig. 37) and should not be confused with the meaning of a super- 
script on a spherical wave function. 
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Since the tangential components of the resulting fields E and H must be con- 
tinuous across the sphere surface, the boundary conditions from equation (13) are applied 
and become, for R = a 


R X g^^^^(rIr’) = RX (r(r’) 

J_ R X V X g1^^^ (R| R’) = Jl R X g1^^^ (r| R') 
^1 ^2 ^ 


Solving equations (23) for the coefficients results in a^^^ = a^^^, b^^^ = b^^^. 

p n ' A n * 

c^^^ = c^^\ and d^^^ = d^^^ where 


^mn ” 


®mn 


Qmn 


.(!)- 


n 


hin(°2> ]' - hinW>]’ 'n<'S> 


Qmn 


(24) 


.(!)- 


n 


hin(^2)I - S h^n<^l>]' jn^^2) 


(25) 


(1) ^nW^ - h^nW)]' ^n<^l> 

[^lhn(^l)]’ in(^2> - [^2W] 


,(D- 


n 


"r [^l^n ^ jn(^2> ' V^l^ 


(26) 


(27) 


Here, 


p-^ = kja 


^2 ~ 
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[pz^(p)] = 1 [pz^ip)] 
= ^2 = "^0 


^1 ■ ^0 



where e^. is the relative dielectric constant and may be complex, and z^p) is either 
or h^(p); also, b^^(p) denotes h^^^(p) since only the spherical Hankel function 
of the first kind appears in this work. 

Substitution of equations (16), (21), and (17) into equations (14) and (15) yields the 
Green's functions for the externally located source. They are as follows: 

Electric source in region I, observation point in region I: 


g^i^)(r|r’) f] y {2- So)ifL±I- (kj) \u' (k. 

^ ^ ^0 ° (n + m) ! \ Jmn [ ®mn ^ 




mn 


+ (kj) 

omn 


qIXih. Qmn 


(28) 


for R > R' , and 


ilr ® 

4t7 

n 

(2 - Sq) 

^ n=l 

o 

II 

a 



(ki)j 

6 

omn j 

omn 


(ki) + 


Ne + (k,) 

rtiTiri 


L 0 


®mn 


(29) 


for a s R < R'. 
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Electric source in region I, observation point in region II: 
ik. 


Gf)(R|R')=!!L f f (2- 


(n - m) ! 
1) (n + m) ! 


M (k^) Mf ) (k ) 


0" 


mn 




Qinn Qinn 


(30) 


for R < a. 


For the source located inside the dielectric sphere, again consider the geometry of 
figure 37, but now with R’ < a. Let R') represent the free- space Green’s 

function in an infinite region with the same constitutive parameters as that of region n. 
From equation (16) 


g1^^^ (R[r’) = ^? V (2 - s„) (n - m)! 

4?7 ' j Z— J " f-k/rs I 1 \ I t 


n=l m=0 


n(n +1) (n + m) ! 


+ (k2) (kg) 


;mn 


;mn 


qIUTI qHiii 


(R > R') (31a) 


gP>(R|R’)=^! 1? y 7^ (n-m)! 

^ 0 n(n + I) (irn^)! 


4 t 7 


n=l m=0 


M (k2)ML^^) (k2) 

Qmn Qmn 


+ N 


(kj) 


mn 


(R < R’) (31b) 


Using the method of scattering superposition, (R| R') is treated as consisting of two 

parts, namely, 


Gp) (rIr’) = G^22) (r|r.) + g {^ 2 ) 

G^^2) ^ 5(12) 


(R^ 

a) 

(32) 

(R ^ 

a) 

(33) 
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where represents the portion of Gg which is scattered from the inside surface 

of the dielectric sphere and is the transmitted field. It is noted that G^^^^ must 

be finite at the origin, that must satisfy the radiation condition at infinity, and that 

the boundary conditions must be met at the sphere surface. Thus 



CO n 


(2- V 


2n + 1 (n - m) ! 
n(n +1) (n + m) ! 


(ki)M’ (k^) 


;mn 


mn 


+ bf ) (k, ) N' (ko) 

Qmn Qinn Qmn 


(R > a) 


and 


S!22),s|5,>,“2 


"3s 


(r|r’)=^ ~ ' c^^^ M (k2>M^ (kg) 

n^l m^ ^ (n + m) ! |_o“” 


„»2> 


Qmn Qmn 


mn 


(R < a) 


(35) 


where the superscript (2) on the scattering coefficients denotes that they pertain to the 
source located in region II. 

The coefficients are now found by applying the boundary conditions at the sphere 
surface, namely, for R = a 

R X (rIr') = R X G^^^^ (r|r') 

J_ R X V X G^^^^ (R| R') = R X V X G^^^^ (r|r') 

Again it is found that a^^^ = a^^\ b^^^ = b^^\ c^^^ = c^^\ and 

qIUFi qIiiii Qtxiri 

for the case 


j(2) _ where 

a ’ 

Qmn 
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( 36 ) 


”~ ■ -■- ■ 

^j(2) ^ [^2jn^^2 j ^n^^2^ ~ hV^^2^] jn^^2^ 

^ [p23n^P2^] ■ sFt [^l*'n^P] 

^(2) _ 

— ^ — . — , — 

h^n^^2^] ■ h'^nW^] ^n^^2^ 

^(2) ^ h^n^l^]* ^1^2) ~ 

hjnM'V^l)- h*^n(^l)lM^2)^r 

Substitution of equations (31), (34), and (35) into (33) yields the desired Green's 
functions for the internally located electric current source as follows: 

Electric source in region n, observation point in region I: 


(R|R’) = !li ^ Y' (2- SQ)jf^-^i.. (kjM’ (k,) 

4. ^ ^0 0 n(n + l){n + m)t « ^ 


Qinn Q mn 


bf n 1'^ ,(kl)N; (kj) 


(R > a) (40) 


Electric source in region n, observation point in region II: 


<x> n 


Z! 0^) 


n=l m=0 


n(n +1) (n + m) ! | [^ Qmn 
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for a ? R > R’ and 


Sp>(R|R')=!^ ^ y (2- wiill /m (k,) (kj) 

1^1 ^0. n(n+l)(n + m)! \ gmn ^ [ «nm ^ 


+ c®> (kj) 

“ Qinn 


+ N. 


Cmn 


0=2) 


'n;« (kj) +dfN' (kj) 


;mn 


mn 


( 42 ) 


for R < R’. 

The dyadic Green's function of the fourth kind may be used to obtain expressions 
for the fields due to the magnetic current source. The magnetic field is given in terms 
of these functions for the appropriate regions of source and observation points as follows: 


H(R) — lo> 


H(R) = iwej 


H(R) = ic^£2 


H(R) = i^ £2 


III 

III 

III 

III 


G^^^^RIR') dV’ 

(R|R») k J^(R') dV 
G^^2^ {B\ R') . Jj^(R') dV 
(R|R’)‘ jj^(R') dV 


(43) 


(44) 


(45) 


(46) 


where is the magnetic current density. 

The functions G^ may be derived from 5g using the principle of duality and are 
as follows: 
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Magnetic source in region I, observation point in region I: 

5f)(R|S')=^ f T (2- <" : /m») (ki) 

■*” 'S ^0 ° n(n + l) (n + m)! 1 gmn 


M; (kj) 


^mn 


+ (kj) 


mn 


id) 


" «a ' “‘l* 

Qinn 


fmn rtinn 


L 0 


( 47 ) 


for R > R’ and 


gP^(rIr') = ^ y] {2 - 8 q ) -HILLL ~ • 


n 


n=l m=0 


n(n + 1) (n + m) ! 


M (kj) 

®mn 


(kj) 

^ nmn 




(kj) + 

Omn 


N. (kj) (kj) 

e — i n 6j^j^ 


:mn 


S;,<>)(ki)l 

omn 


for a < R < R'. 


Magnetic source in region I, observation point in region II: 


(48) 


gf)(R|R.)=i^ f f (2- 


(n - m)t 
1) (n + m)! 


td^>M (ko)M;^"^ (kj) 
^ ®mn ®mn ^ 




y^ixin. /%mn 


(49) 


for R g a. 

Magnetic source in region II, observation point in region I: 


ikn 


=-i 2^ 7] ■ 

n=l m=0 


Sn) 


2n + 1 (n - m) ! 
n(n + 1) (n + m) 


-R 


(kj)M' (k„) 
®mn ^ ®mn ^ 


for R > a. 
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n e 
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(k.) 


mn 




(kg) 


mn 


(50) 
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Magnetic source In region II, observation point in region II: 


Gf2)(R|R.) y (2 - Sq) ^ K (k,) 

^ fcl ifel n(n + 1) (n + m)! Qinn ^ 


.d(2) 


Qinn J 


qIHII ^ qIhii qHIii qUih j 


for a > R > R', and 


Sf >(r|r') = y y ) 2^ (n - m)! (m , [«.(!) (k 

fcl ^0 ^ n(n + l)(n + m)! ^mn |_ ^mn 


0 


+ N. 


;mn 


(kj) 


«£’ <‘2) " 'n ’ N' (kj) 

/-\tnn /-iTXin 


LO 


(52) 


for R < R'. 


Electric Fields Due to a Huygens' Source in the Presence 
of a Dielectric Sphere 

For the source located as shown in figure 36, the electric and magnetic current 
densities are represented by 


5(R’ - bj) 8(^') S(<^’) . 

J(R') = C„ i X (53) 

bj sin 

and , 

_ _ s(R' - b.) s(e') sOT 

•'«»’> = - ^ y (54) 

bj sin d' 

respectively. The electric fields due to the current given in equation (53) are easily found 
by using equations (8) to (11). In order to obtain the electric fields due to the magnetic 
source, however, equations (43) to (46) are used together with Maxwell’s equation for a 
source-free region 
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(55) 


E =-l VX H 

The curl operator is easily performed in view of the symmetrical relation for the vector 
wave functions of equation (20). 

The resulting fields for the combined source may be represented as follows; 
Source in region I, observation point in region I: 




for R > bj^ and 


Eh(R) = 






^ 2n + 1 
~~i n(n + 1 


4 t 7 - ^ 




, hV'-s)]') 

- 1-t 


{“Olnfl'l' 


(56) 




for a < R < bj^. 


Source in region I, observation point in region II: 


for 0 I R ^ a. 


E {{W - “oinO^l 

n=l ^ ' 


On ’ Neln(^) 
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Source in region n, observation point in region I: 


n=l ^ ^ 


Source in region II, observation point in region II: 


1h(r) = :^ y 


r'^^e 2n+l J. , X 

4 ^ ^ + 1) \V 


- i c^iMlrM(i) 


+ + 


hin(^4)]’ 




for as R > bj. 


Eg(R).:^ f ^ 
^ 4tt n(n + 1) 
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Nn<^ 4 )]’ 
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ih„(P 4 )+iC V^ 4 )rNeln<k 2 ) 


for R < bj, where 


A?’ ' i„h) " 4‘* W 


bin<P3>]' ^ 1,(1) b^PS*]' 


o(l) = ^ jj(l) 

n n„ n 



(64) 


p(l) _ ,,{1) 

n '■ ,3 n 

P3 = kjbj 

A (2) - i (n ) [T~ 


(65) 


(66) 


(2) _ h^rS^^] . (2) 

n - 

(67) 


(68) 

,(2)_i ['’4l»('’4)]' ,(2)^ 

(69) 


% = k2b^ 

The coefficients in equations (56) and (57) were combined into and 

D since the electric field far from the antenna is of primary interest. Here, A and 
n 

represent the field contribution from the electric current source and, similarly, 
and were derived from the magnetic current source. Again, jthe superscript on the 
coefficient denotes the source region. 

For the source located on the sphere surface, these coefficients may be reduced 
using the relation 

Jn(^) = ? 
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Hence 
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= A^2) 


a(i)=A 
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’ 
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(70) 


(71) 


(72) 


(73) 


In the far zone the spherical Hankel functions of the first kind may be replaced by 
the leading term of their asymptotic expressions. Thus 


h„(kR) - (-i)'"''^ ® 


ikR 
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(74) 


For large values of kR the spherical wave functions become 


(kR) - (-i)"^^ - 

gmn kR Qmn 
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kR Qmn 
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where 


and 
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;mn 
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Using these results in equations (56) and (59), the far-zone field expressions are 
obtained for the Huygens' source in the presence of a dielectric sphere 


where 
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(79) 


(80) 

(81) 


(82) 
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for the source located in region I or II, respectively. 

In order to recover the far-zone field expressions due to the electric or magnetic 
source only, ~ ^n “ ^ ^n ~ ®n “ ^ equations (79) to (83), 
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Properties of a Huygens' Source in the Presence of a Dielectric Sphere 

Radiated power .- An expression for the total radiated power W is now derived by 
integrating the total power flow across an infinitely large sphere, centered at the origin, 
using the relation 

W = — f f E • E* dn (84) 

2^0 •^0 -^0 

where dn = sin d de d<p and the asterisk denotes complex conjugation. 

The various scalar products formed in equation (84) are now investigated. It is seen 
from equations (77) and (78) that each term of the product contain a factor 

of the form 


pJ^(cos $) aP^ (cos e) 


(85) 


sin0 


W 


Using the differential relations for the associated Legendre functions (ref. 5), it is found that 


P^(cos &) aP^ (cos 0) 2 


sin & 
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It now follows that 
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sin 9 d 0 = 0 
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since 
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P™ (cos 0) P^(cos 9) d0 =0 
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Thus 
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(89) 


Also, performing the integration results in 


Jj"'oin-“ou =JKln-"eM dn = 

2n + 1 
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(nj^i) 

(n = 1) 

and where the formula (Stratton, ref. 5) 
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2n + 1 (n - m) ! 




(92) 


has been used. Inserting equation (79) into equation (84) and evaluating the integrals gives 
the desired relation for the radiated power as follows: 


W = 


n 

2 




CD 



n=l 


(93) 


where Vq = 1207t ohms. 


Radiation resistance .- Suppose that instead of having an infinitesimal point source, 
the electric current element were of a short lei^h i with a constant magnitude Iq. 
Then = Iq^, and the radiation resistance for this system is easily obtained. 

The ra.diation Qn an+onn<i to oq 

“ w «h^4.AW«^ J.u.xc^ AK., V*X^A***^WC 4%D 


Rr = 


2W 


(94) 
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That is, assuming no losses, a current element with an rms value of Iq flowing on an 
antenna with a radiation resistance R^., will radiate an average power W. In order to 
determine the radiation resistance of the source in the presence of the dielectric sphere, 
equation (93) is substituted into equation (94). Hence 


= 15(k^^)2 (2n + 1) {a^a^ + /3^*) 

n=l 


(95) 


Directivity .- The directivity in the forward direction relative to an isotropic source 
is defined by 


4t 7 (Radiation intensity in g = 180^ direction) 
Total power radiated 


(96) 


since for the sphere concentric with the origin, and the source located on the +z-axis, 
the relative field intensity in the -z {e = 180*^) direction is of primary interest. 

From equations (77) and (78), it is evident that 
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Using this result in equation (94) an expression for the radiated power intensity in the 
-z -direction is obtained. Thus 




E ♦ E* 


15K^C^ 
I e 
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^(i)Nan./3n) 
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Combining equations (101), (93), and (96) the desired directivity formula results; 
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□0 

^ ' 2ri[ + 1 /*\n / \ 
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( 102 ) 


2n + 1 




For comparison, it will be of interest to compute the directivity of a uniformly 
illuminated circular aperture. The normalized field pattern of a large uniform aperture 
can be found by the application of Huygens' principle and is given by Kraus (ref. 6) as 


/ 77-d 

OiJi \ \ sin^ 

E(^) =1^ J_AA^_ 

77 d sin 8 


(103) 


where d = 2a is the aperture diameter. For this case, the directivity is given by 


n - IE(0)I^ 

^0 



and hence, 



Jq sin 8 
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Efficiency .- For a lossy dielectric, is complex and may be represented by 


^2 = e' + U" 


(106) 


where e' and e" are the real and imaginary parts of the permittivity, respectively. 
The conductivity a is related to e” by 


cr= 


(107) 


and the loss tangent is defined as 

e” 

tan S = -- 
e’ 

Equation (106), therefore, may be cast in the form 
= ^'(1 i tan s) 


(108) 


(109) 


The propagation constant k 2 will also be complex, being equal to 
k 2 = k’ /I + i tanT 


( 110 ) 


Thus in the previously derived equations « 2 ^2 simply replaced by their com- 

plex values to account for the dielectric loss. 

To find the power dissipated within the sphere, the source is placed in region I close 
to the surface of the sphere. With the electric field within the sphere may 

be expressed by 


Ejj(R) = 


-kj^ai^C^ 


E 


2n + 1 


in ^ n(n + 1) 




M 


01n^2^ +^'^eln 


(kg)] 


( 111 ) 


where 




PVcos 0) 


n 


cos<^ - 


3P^(C0S 0 ) 


Sin <M>\ 


sin (9 


d6 


( 112 ) 
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^ n(njJO P^(cos 9) cos 0R 


1 r. i' ^Pn^COS( 9 ) ^ P'[(cos 0) 

+ — Loj_(p) cos(?b0 sin 

^ ^ L Be sine 


P = k2R 




\/V ^3 


= d(«LlA_LL.iJ!_h„(P3) 


and and dj|^^ are given by equations (26) and (27), respectively. 

The power dissipated within a region of conductivity a- is given by the formula 


’D'il/J 


crE • E* dV 


For a homogeneous sphere of radius a this equation becomes 


,2-rr -n a 


Pd = 2 J / J 


0 0 0 


Because of the orthogonality of the associated Legendre functions and their 
derivatives. 


f ’^Oln * ^el^ d5 - 0 
•"O 
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Also, it is easy to show that 



^01 ' ^5 It 


27rn^ 


(n + 1)^ 
2n + 1 


i„Mi*W *n{ 


(117) 


and 





_ 2770^ (n + 1)^ 
P^*(2n -I- 1) 



+ 


[pinW] 


where equation (92) and the formula (also given by Stratton in ref. 5) 


(118) 


7T 

f* 

(cos 0 ) (cos 0 ) sin 0 de 

0 


2 (n + m) ! j, 

2n + i (n - m)l ^ 


(119) 


have been used. With the aid of equation (118) the 0 and 4 > integrations are performed 
and an expression for the dielectric loss is obtained. Thus 


Pd = 


,2 2 2^2 

OK^OJ 

16t7 


^ (2n + 1) r l)in^p)3n((°) 

n=l '’o I PP ^ 


+ 





dR 


(120) 


or 




2|'rl 


CP ^ 

J](2n + 1) r |vpj>)|^+n(n + l)|dJp)l2 + t.^[pj^(p)]'|2dR 
n=l h 


( 121 ) 


Now that expressions for the dissipated power Pjj (eqs. (120) and (121) and the 
radiated power W (eq. (93)) are available, the antenna efficiency may be computed from 
the relation 
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Radiated power W 

= (122 

Total input power W + Pj^ 

NUMERICAL RESULTS 

The far-zone electric fields radiated from a Huygens’ source on the surface of a 
dielectric sphere (see eq. (79)) were numerically calculated using the computer program 
listed in the appendix. Since the series for the far field converges in n in a manner 
similar to the spherical Bessel functions, it was found that N § kD + 7 terms provided 
sufficient accuracy in the calculations within plotting error. 

The normalized radiation patterns for various diameter spheres, computed in 
1*^ increments, are presented in figure 38, Here the E- and H-planes correspond to 
the <^1 = 0 and 4> = tt/ 2 planes, respectively. In each case a dielectric constant of 2.57 
with a loss tangent of 0.0065 was used, which is typical of Plexiglas in the frequency 
range of interest (ref. 7). It should be noted that with the source oriented as shown in 
figiu-e 36, the maximum radiation occurs in the -z-direction, and hence the top of the 
computed patterns correspond to 0 = 180^. 

In order to investigate the effect of moving the source away from the sphere, the 
curves of figure 39 were calculated using equation (79) in 2® intervals. Since bj is 
the distance from the origin to the source and a is the sphere radius, the ratios of 
bj^/a shown in figures 39(a) to (h) correspond to the source displaced 0.0, 0.25, 0.5, 
and 1.0 free-space wavelengths from the sphere surface. In a similar manner, the 
radiation patterns of figure 40 correspond to the source spaced 0.1 and 0.5 inside the 
sphere surface. 

By increasing the dielectric constant, directive patterns can be produced for the 
source lying inside the dielectric sphere. Typical patterns computed for 12.70-cm 
spheres operating at 10 GHz are shown in figure 41. 

Figure 42 shows the optimum ratio of soimce to sphere radius for various diameter 
spheres. These points were obtained by fixing the dielectric constant at 4, 6, 9, and 12, 
and varying the source position until the maximum directivity was found. The additional 
points corresponding to the source on the sphere surface were obtained by adjusting the 
dielectric constant. In each case a loss tai^ent of 0.0065 was used. For comparison, 
the upper curve shows the location of the geometrical optics paraxial focus. 

Figure 43 shows the directivity (computed from eq. (102)) of the Huygens' source on 
the surface of the sphere. Again, a dielectric constant of 2.57 and loss tangent of 0,0065 
have been used. For comparison, the directivity of a uniformly illuminated circular aper- 
ture (eq. (105)) is also shown. The abscissa is measured in aperture diameter or sphere 
diameter in free-space wavelengths. The circles and squares in the figure represent 
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directivities derived from measured patterns for waveguide excited Plexiglas spheres 
using the empirical formula 


= -1L (123) 

where K is a constant and 8 ^ and 8 ^ are the half-power beamwidths in the two principal 
planes measured in degrees. By comparing the computed directivity with the beamwidths 
of the Huygens’ source patterns it was found that K = 26 500 yielded an error less than 
0.5 dB for sphere diameters in the range of interest. Since the measured and computed 

patterns have the same general shape, this value of K was used in the above formula. 

The roughness in the curve is caused by a resonance effect within the sphere which has 
been damped by the dielectric loss. Figure 44 was constructed by drawing smooth curves 
through maximum points on directivity curves for other dielectric constants. 

The antenna efficiency, computed using equation (122), is plotted in figure 45. For 
this computation, the source was placed O.OUQfrom the sphere surface to avoid the 
problem of infinite loss associated with a point source in contact with a lossy medium 
(Tai, ref. 8). The dips in the curve at the larger diameter spheres are, again, due to the 
damped resonances. 

Figure 46 shows the radiation resistance, computed from equation (95), together 
with reflection coefficient curves measured for 5.08- and 7.62 -cm Plexiglas spheres. 

It is noted that both the computed and measured curves tend to have maximums spaced 
about 0.25 apart. 


DISCUSSION 

Computed results and measurements using a circular waveguide feed are presented 
in figure 47 for a variety of Plexiglas sphere diameters. These results provide a direct 
comparison of the theoretical Huygens’ source model and the experimental measurements. 

It can be seen by inspecting figure 47 that the agreement between results is good, particu- 
larly down to the -10-dB pattern level. The theoretical patterns do, however, show sharper 
nulls in the main beam and stronger side lobes in the rear direction than do the measured 
patterns. These differences may be due to the scattering from the feed waveguide structure 
used in the measurements or due to the fact that the Huygens' source model does not have 
the same forward directivity as the open-end waveguide feed antenna. Also precise agree- 
ment may not be possible in particular cases because of the lack of knowledge as to where 
to place the position of the equivalent Huygens' source. For example, in the case of the 
AED antennas using spheres, part of the dielectric protrudes into the feed waveguide, 
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particularly for a circular waveguide feed using small spheres. The effects of moving 
the Huygens* source inside and outside the sphere are given in figures 39 and 40. Although 
the patterns are not given here, it has been found (ref. 3) that extending the Huygens* point 
source into curved- line sources, subtending about the same at^le as the waveguide aperture, 
did not significantly alter the calculated patterns. 

As noted above, the radiation patterns tend to degrade as the source is moved away 
from the Plexiglas sphere, A similar result was also seen in the computed patterns, for 
the Plexiglas dielectric constant, as the source is moved into the sphere. These results 
indicate, therefore, that the optimum source location for Plexiglas is near the sphere 
surface. As seen from the curves given in figure 42 and the patterns given in figure 41, 
it should be possible to produce directive patterns with a source located interior to the 
sphere by increasing the dielectric constant appropriately. 

The directivity curve, figure 43, computed using the Huygens* source model, shows 
close agreement with the directivity computed from experimental measurements, and 
confirms that the directivity of the AED antenna is indeed greater than that of a uniformly 
illuminated aperture of the same cross-sectional area. However, gain measurements, 
which include possible losses in the dielectric material, show that such losses are sig- 
nificant, figure 48. Such a behavior can be inferred from the radiation-efficiency curve 
given in figure 45, which clearly shows a loss in gain with increased volume of dielectric 
material, about the order of the measured values in figure 48. 

The computed antenna efficiency decreases with increased sphere size because the 
propagation path through the lossy dielectric material is increased. It is noted, however, 
that near certain frequencies the efficiency curve exhibits large changes in efficiency. 
Similar chaises, although not as pronounced, occur in the directivity curves. These 
changes are caused by a resonance effect within the dielectric sphere, such resonances 
being damped by losses in the material. For example, consider the computed radiation 
patterns given in figure 49. For the lossless case at resonance, much of the energy is radi- 
ated into side and back lobes, resulting in a drop in directivity and changes in the radiation 
resistance. With increased loss the side and back lobes are reduced at a sacrifice in 
efficiency. This resonance effect becomes more pronounced with increased sphere size 
and dielectric constant. Further detailed discussion of the resonance effect is given in 
reference 3. 


CONCLUDING REMARKS 

The results of a study of the radiation from waveguide- excited uieleetric objects 
have been presented. The radiation patterns measured for small Plexiglas cubes and 
spheres exhibit gains in excess of those which can be obtained using an optimized horn of 


40 



the same cross-sectional area. It was also noted that the depolarization with the dielectric 
loading was no greater than that of an open-end waveguide. 

By using the model of a waveguide opening onto an infinite dielectric slab-covered 
ground plane, it was found that the input admittance of the aperture excited dielectric 
(AED) antenna could be roughly predicted. 

The coupling between adjacent loaded antennas was reduced below that of the unloaded 
ivaveguide aperture. This was in agreement with near-field amplitude measurements taken 
around the outside of a waveguide- excited sphere, which showed that most of the energy 
was concentrated in the forward part of the dielectric object. 

Expressions for the radiated fields and other antenna parameters were derived for 
the Huygens' soxirce in the presence of a dielectric sphere, using the method of dyadic 
Green's functions. The calculated radiation patterns and directivity using this source 
model on the surface of the sphere closely corresponded to the measurements. A dielec- 
tric constant of about 3.0 tended to yield the maximum directivity for spheres greater 
than 2.5 wavelengths in diameter. It has been confirmed from the methods of geometrical 
optics (ref. 3) that this value of dielectric constant yields little phase error across an 
equivalent aperture plane. Desirable patterns could also be produced, with the source 
located inside the sphere, if the dielectric constant were increased appropriately. 

Finally, it was noted that the resonance effect of the dielectric sphere became 
stronger with increasing dielectric constant and sphere diameter. With dielectric materials 
such as Plexiglas, however, the resonance does not degrade the radiation pattern, but is 
manifest as a decrease in the antenna efficiency. This is due to the damping effect of the 
dielectric loss. 

Lar^ley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 21, 1973. 
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APPENDIX 


COMPUTER PROGRAM 

Although numerous computer programs were used throughout this research, the one 
of primary significance, together with the necessary subroutines, is listed in this appendix. 
The program was written in standard IBM FORTRAN IV for use with the G level compiler 
and was run on an IBM System 360/67 computer. 

The MAIN program calculates and prints out the normalized far-zone radiation 
patterns (magnitude, phase, and intensity in dB) in the two principal planes for a Huygens’ 
source on the surface or inside a lossy dielectric sphere. It also prints the coefficients 
of the functions m and n and the directivity and provides the pattern data (truncated 
to -40 dB) in an array suitable for plotting. This program may easily be modified to com- 
pute data for only the electric or magnetic current source or, for the case where the 
source lies outside the dielectric sphere, to make the appropriate changes in the coeffi- 
cient formulas of lines 35 to 39. 

The input parameters to the MAIN program are defined as follows: 

FORTRAN parameter Description 


DIAM Sphere diameter, \q 

PER Dielectric constant 

TAN Dielectric loss tangent, XIOOO 

SP Distance from source to sphere surface, 

thing Increment in e variable, degrees, (The parameter THING 

should divide evenly into 180°. If THING = 1 , for exam- 
ple, data are calculated from 0 = 0° to 0 = 180° in 
1° intervals.) 

Subroutine ASLEG is used to generate arrays of associated Legendre functions of the 
first order and their derivatives. The parameters are as follows: 
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APPENDIX — Continued 


FORTRAN parameter Description 

NMAX Maximum value of degree, N, desired 

THING Increment in d variable, degrees. (See note above.) 

Names of single-precision two-dimensional arrays in which 
the output functions are returned 

pJj(cos M) 

where M = (N - 1)THINC 

sin M 

BP^(cos M) 

— , where M = (N - 1)THINC 

Subroutine CBESS generates arrays of spherical Bessel and Hankel functions of 
complex arguments and their derivatives. The parameters are as follows: 

FORTRAN parameter Description 


SLEGEN 

dlegenJ 

SLEGEN(N,M) 

DLEGEN(N,M) 


BES 

DBES 

HANK 

DHANK 


Names of single-precision complex arrays in which the out- 
put functions are returned 


BES(N) 

Jn-i® 

DBES(N) 


HANK(N) 


DHANK(N) 


X 

Double precision complex argument 

NMAX 

Maximum value of order (N) desired 


Double precision arithmetic is used throughout these subroutines, with the final out- 
put converted to single precision. 
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TABLE 1.- MEASURED GAIN OF AED ANTENNAS USING CIRCULAR 
OR SQUARE WAVEGUIDE FEED 


Size and shape of 
dielectric object 

Gain, dB, using circular waveguide 
measured at — 

Gain, dB, using 
square waveguide 
at f = 10.0 GHz 

5.4 GHz 

5.5 GHz 

5.6 GHz 

Open-end waveguide 

9.45 

9.12 

8.33 

10.00 

5. 08 -cm cube 

12.70 

12.95 

12.54 

13.82 

5.08-cm sphere 

10.32 

10.75 

9.92 

13.36 

7.62-cm sphere 

13.05 

13.20 

13.00 

17,40 

10.16-cm sphere 

17.26 

17.02 

16.45 

19.72 


TABLE 2.- MEASURED VSWR OF A CIRCULAR WAVEGUIDE, 

5. 08- cm SPHERE 


Frequency, 

GHz 

Open-end 

waveguide 

Without 

taper 

VSWR for - 

Taper (a) 

Taper (b) 

Taper (c) 

5.0 

1.95 

2.40 

1.26 

1.68 

1.70 

5.2 

1.55 

2.50 

1.68 

1.08 

1.17 

5.4 

1.08 

2.40 

1.32 

1.40 

1.45 

5.6 

1.20 

2.42 

1.13 

1.08 

1.01 

5.8 

1.13 

2.40 

1.16 

1.02 

1.21 

6.0 

1.52 

2.20 

1.63 

1.50 

1.53 
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TABLE 3.- COUPLING OF PLEXIGLAS- LOADED CIRCULAR WAVEGUIDE ANTENNAS 

AS A FUNCTION OF FREQUENCY 

























Figure 2.- Photograph of the circular waveguide feed with Plexiglas objects. 
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(d) 7.62-cm sphere. (e) 12.70-cm sphere. 

Figure 3.- Radiation patterns of AED antennas in the E-piane for a circular waveguide feed at 5.0 GHz. 




(d) 7.62-cm sphere. (e) 12.70-cm sphere. 

Figure 4.- Radiation patterns of AED antennas in the H-plane for a circular waveguide feed at 5.0 GHz. 






(a) Open slot, 


{b) 5.08-cm cube 


(c) 5.08-cm sphere 


(d) 7.62-cm sphere. 


(e) 12.70-cm sphere. 


Figure 6.- Radiation pattern for AED antennas in the H-plane for a slot feed at 10.0 GHz. 





(d) 7.62-r,m sphere. (e) 12.70-cm sphere. 

Figure 7.- Radiation patterns for AED antennas in the E -plane for a rectangular wave^ide feed at 10.0 GHz 
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(a) Open end. 


(b) 5.08-cm cube. 


(c) 5,08- cm sphere. 
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(d) 7.62-cm sphere. (e) 12.70-cm sphere. (f) 15,24-cm sphere. 

Figure 9.- Radiation patterns for AED antennas in the E-plane for a square-horn waveguide feed at 10.0 GHz 
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(e) 12.70-cm sphere. 











(d) 7.62-cm sphere. (e) 12.70-cm sphere. (f) 15.24-cm sphere. 

Figure 10.- Radiation patterns for AED antennas in the H-plane for a square-horn waveguide feed at 10.0 GHz 






(a) E-plane; 0.0-cm spacing. 


(b) E-plane; 0.63 5-cm spacing. 


(c) E-plane; 1.27-cm spacing 
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(d) H-plane; 0.0-cm spacing. (e) H-plane; 0.635-cm spacing. (f) H-plane; 1.27-cm spacing. 

Figure 11.- Radiation patterns of a circular waveguide feed exciting a 12.70-cm Plexiglas sphere 
placed various distances from the waveguide aperture operating at 5.0 GHz. 
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Figure 12.- Photographs of hemispherically capped and rectangular block AED antennas. 


(b) Rectangular. 
Figure 12.- Concluded. 


L-73-3799 



(a) E-plane; 5.08 x 5.08 x 2.54 cm. (b) E-plane; 5.08 x 5.08 x 5.08 cm. 



(c) H-plane; 5.08 x 5.08 x 2.54 cm. (d) H-plane; 5.08 x 5.08 x 5.08 cm. 

Figure 13.- Radiation patterns of a circular waveguide exciting a 
Plexiglas block of different dimensions at 5.0 GHz. 
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(d) H-plane; 0.0 cm. (e) H-plane; 1.27 cm. (f) H-plane; 2.54 cm 


Figure 14.- Radiation patterns of a circular waveguide antenna exciting a 5.08-cm hemispherically 
capped circular Plexiglas cylinder at 5.0 GHz as a function of cylinder length. 
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(b) Photograph. 
Figure 15.- Concluded 



(a) E-plane; inverted position. (b) E-plane; direct position. 



(c) H-plane; inverted position. (d) H-plane; direct position. 


Figure 16.- Radiation patterns of a circular waveguide exciting a shaped 
Plexiglas obstacle at an operating frequency of 5.0 GHz. (The dashed 
pattern is for the unshaped sphere in the same waveguide.) 
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Beamvidth, deg 









Frequency, GHz 


Figure 18.- Gain of spherical AED antennas fed by a circular waveguide compared 
with the gain of the optimum horn having the same cross section. 











Theory 



Frequency, GHz 

Figure 21.- Calculated and measured values of the VSWR of 
circular waveguide fed AED antennas over the frequency 
range from 5.0 to 6.0 GHz. 
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(a) Taper (a). (b) Taper (b). (c) Taper (c). 

Figure 22.- Sketch and photograph of tapered Plexiglas sections for matching the 

circular waveguide AED antennas. 
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(d) Photograph of tapered Plexiglas sections. 
Figure 22.- Concluded. 
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Figure 23.- Sketch showing the detection circuit for the near-field amplitude 
measurements. {Diode, MA 41512; Zj^, 1000 ohms; Z2, 1000 ohms; 

Z3, 500 pF; Ri, 10 000 ohms.) 
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Figure 24.- Geometry of the near-field amplitude measurements using 
a dipole probe and a circular waveguide feed. 
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Figure 25.- Continued. 
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(a) E-plane at Y = 0.508 cm, X = 0.0 cm, and f = 6.0 GHz. 

Figure 26.- Near-field amplitude measurement of a circular waveguide feeding a 5.08-cm Plexiglas cube. 



(b) H-plane at Y = 0.508 cm, Z = 0.0 cm, 
Figure 26.- Continued. 
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Y,cm 

(c) Along the side at X = 3.048 cm, Z = 0.0 cm, and f = 6.0 GHz, 

Figure 26,- Concluded. 




Z, cm 

(a) E-plane at Y = 0.508 cm, X = 0.0 cm, and £ = 6.0 GHz. 

Figure 27.- Near-field amplitude measurement of a circular waveguide feeding a 5.08-cm Plexiglas sphere 
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(b) H-plane at Y = 0.508 cm, Z = 0.0 cm, and f = 6.0 GHz. 
Figure 27,- Continued, 
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(a) E-plane at Y = 0.508 cm, X = 0.0 cm, and f = 6.0 GHz. 

Figure 28.- Near-field amplitude measurement of a circular waveguide feeding a 7.62-cm Plexiglas sphere. 







(a) E -plane at Y = 0.508 


Figure 29.- Near-field amplitude measurement of a < 
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Z = 0.0 cm, and f = 6.0 GHz. 

jular waveguide feeding a 10.16-cm Plexiglas sphere. 
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(a) E-plane at Y = 0.508 cm, Z = 0.0 cm, and f = 6.0 GHz, 


Figure 31,- Near-field amplitude measurement of a circular waveguide feeding a 15.24-cm Plexiglas sphere. 
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(c) Along the side at X = 8.128 cm, Z = 0.0 cm, and f = 6.0 GHz. 

Figure 31.- Concluded. 
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(a) E-plane: 0.95^. 


(b) E-plane; 2.22k, 


(c) E-plane; 6.03X 






(d) H-plane; 


(e) H-plane; 2.22^. 


(f) H-plane; 6.03^. 


Figure 32.- Radiation pattern of a circular waveguide antenna radiating in the presence of an identical antenna 
operating at a frequency of 5.0 GHz and at various center-to-center spacings. (The active antenna is on 
the 270° side of the pattern.) 







(a) E-plane; 0.95\. 


(b) E-plane; 2.22k. 


(c) E-plane; 6. 03 A. 





(d) H-plane; 0.95X. 


(e) H-plane; 2,22x. 


(f) H-plane; 6,03X. 


Figure 33,- Radiation pattern of a circular waveguide antenna loaded with a 5.08-cm Plexiglas sphere radiating 
in the presence of an identical loaded antenna, operating at f = 5.0 GHz, and at several center-to- center 
spacings. (The active antenna is on the 270*^ side for 0.95X and 2.22,v and on the 90° side for 6.03x.) 
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Figure 34.- Photograph of a pair of circular waveguides located in a ground plane which 

were used for coupling measurements. 
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Coupling in dB down Coupling in dB down 


^ Open-end. waveguide 
I ^.06-cm cube 



5.08 7.62 10,16 12.70 15.24 17.70 

Cent er-to-c enter spacing between waveguides, cm 



Center-to-center spacing between waveguides, cm 

Figure 35.- Coupling between a pair of C-band circular waveguides loaded with 
dielectric objects for two alinements of the electric field; f = 5.8 GHz. 
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Figure 36.- Geometry for a Huygens' source in 
the presence of a dielectric sphere. 


Observation point 



Figure 37.- Geometry for a current source in the presence 
of a dielectric sphere. 
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(a) 7. 62 -cm sphere at 

5.0 GHz; E -plane; 

D = 1.27X. 



(b) 7, 62 -cm sphere at 

5.0 GHz; H-plane; 

D = 1.27>^. 



(c) 12.7-cm sphere at 

5.0 GHz; E -plane; 
D = 2.12\. 


(d) 12.7-cm sphere at 

5.0 GHz; H-plane; 
D = 2.12V. 


Figure 38.- Calculated radiation patterns using a Huygens' source 
on the surface of a Plexiglas sphere. 





(a) E-plane; bj/a = 1.0. (b) H-plane; bj/a = 1.0. 



(c) E-plane; bi/a=1.15. (d) H-plane; bj/a = 1.15, 

Figure 39,- Radiation patterns for a source located outside the sphere. 
D = 3.39X; a 2.57; tan S = 0.0065. 
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Figure 42.- Source location as a function of dielectric constant 
for optimum directivity. 
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^Uniform aperture 


Measured pattern (circular waveguide) 
Measured pattern (rectangular waveguide) 
Computed from equation (102) 
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Diameter, wavelengths 


Figure 43.- Directivity of the Plexiglas sphere. Sj. = 2.57: tan g = 0.0065. 
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Figure 44.- Directivity of spheres with tan S = 0.0065 for a 
range of dielectric constants. 
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Figure 46.- Calculated radiation resistance and measured swept frequency 
reflection coefficients for Plexiglas spheres. 
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(a) D = 1.27\. 

Figure 47.- Comparison of calculated and measured radiation patterns 
for a circular waveguide antenna exciting a Plexiglas sphere of 
various diameters. 












